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ABSTRACT 



Aims. We present a study of the diffuse X-ray emission in the halo and the disc of the starburst galaxy NGC 253. 

Methods. After removing point-like sources, we analysed XMM-Newton images, hardness ratio maps and spectra from several regions in the 
halo and the disc. We introduce a method to produce vignetting corrected images from the EPIC pn data, and we developed a procedure that 
allows a correct background treatment for low surface brightness spectra, using a local background, together with closed filter observations. 
Results. Most of the emission from the halo is at energies below 1 keV. In the disc, also emission at higher energies is present. The extent of 
the diffuse emission along the major axis of the disc is 13.6 kpc. The halo resembles a horn structure and reaches out to ~9 kpc perpendicular 
to the disc. Disc regions that cover star forming regions, like spiral arms, show harder spectra than regions with lower star forming activity. 
Models for spectral fits of the disc regions need at least three components: two thermal plasmas with solar abundances plus a power law and 
galactic foreground absorption. Temperatures are between 0.1 and 0.3 keV and between 0.3 and 0.9 keV for the soft and the hard component, 
respectively. The power law component may indicate an unresolved contribution from X-ray binaries in the disc. The halo emission is not 
uniform, neither spatially nor spectrally. The southeastern halo is softer than the northwestern halo. To model the spectra in the halo, we needed 
two thermal plasmas with solar abundances plus galactic foreground absorption. Temperatures are around 0.1 and 0.3 keV. A comparison 
between X-ray and UV emission shows that both originate from the same regions. The UV emission is more extended in the southeastern halo, 
where it seems to form a shell around the X-ray emission. 

Key words. X-rays: galaxies - X-rays: ISM - Galaxies: individual: NGC 253 - Galaxies: halos - Galaxies: ISM - Galaxies: starburst 



1. Introduction 

The diffuse X-ray emission of starburst galaxies can be quite 
prominent. Especially in galaxies that we see edge-on, we 
can find very complex emission from galactic halos. One fa- 
mous example is the starburst galaxy NGC 253 in the Sculptor 



Group . It is close enough (2.58 Mpc, l'=750 pc. iPuche et al 



199 II) to resolve structures in the disc and halo, and to sep- 
arate the detected point-like sources from the diffuse emis- 
sion. Also, it is seen almost edge-on (78.5°, |Pencdll980l) . 
so an unobscured analysis of the halo emission is possible. 
NGC 253 has been observed in X-rays many times. Th ere are 
observations with Einstein ( e .g. iFabbiano & Trinchieri 1984). 
ROSAT (e.g. IPietschl 1 19921: iRead et all 1 19971: iDahlem et al 
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19981: IVogler & Pietschl 1 19991: 



(e.g. iPtak et al 



Pietsch etail l2000h 

ffl, BeppoSAX (e.g . ICappi et al 

12001b iBauer et al 



ASCA 



1999) 



2007) 



2002, 



XMM-Newton (e. g. IPietsch et al. .. 
and Ch andra (e.g. IWeaver et al.ll2002t IStrickland et al. 
2004allbl) . While with some instruments one was not able 
to separate emission from point sources and diffuse emis- 
sion, other instruments, especially ROSAT, XMM-Newton, and 
Chandra, do have a narrow enough point spread function to do 
so. We here report on the first extensive analysis of the diffuse 
emission in NGC 253 with XMM-Newton. 



2. Observations and data reduction 



NGC 253 was observed with XMM-Newton dJansenet al.l200ll) 
during three orbits in June 2000 and June 2003, using all 
of the European Photon Imaging Camer a (EPIC) instruments 
dStriider et alj 1200 U iTurner et al.l 12001 ), the two co-aligne d 
RGS spectrometers, RGS1 and RGS2 dden Herder et"ai1l200lh . 
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Table 1. XMM-Newton EPIC pn NGC 253 observation log. 



Nr. 


Revolution 


Obs. ID. 


Obs. Dates 


Pointing 


Direction 


PA. 


Filter 


T 

L exp 


-Lexp, clean 




Number 






RA/DEC (12000) 


(deg) 




(ks) 


(ks) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


89 


0125960101 


2000-06-03 


00:47:36.74 


-25:17:49.2 


56.9 


Medium 


60.8 


24.3 


2 


89 


0125960201 


2000-06-04 


00:47:36.57 


-25:17:48.7 


57.0 


Thin 


17.5 


3.1 


3 


186 


0110900101 


2000-12-14 


00:47:30.20 


-25:15:53.2 


233.8 


Thin 


24.4 


4.4 


4 


646 


0152020101 


2003-06-19 


00:47:36.89 


-25:17:57.3 


53.8 


Thin 


113.0 


47.9 



and the Optical Monitor (IMason et al.l 1200 II) . for a total of 
about 216 ks. The revolution number, observation identifier, 
observing date, pointings and orientation of the satellite (P. A.), 
and the total exposure times for the EPIC pn camera (T exp ) are 
shown in Table Q] 

Throughout the following analysis, we used the MOS data 
only to detect and remove point-like sources. We were espe- 
cially interested in low-surface brightness diffuse emission at 
energies below 1 keV, where the MOS detectors have a lower 
sensitivity than the EPIC pn. By not utilising the MOS data for 
the analysis of the diffuse emission, we avoided a higher back- 
ground noise level. 

We analysed the data using the Science Analysis System 
(SAS), version 7.0.0. In a first step, we cleaned the EPIC pn and 
MOS observations by excluding times, where the count rate 
over the whole detector above 10 keV exceeded the count rate 
during quiescent times. This cleaning was done to avoid times 
with high particle background, which result in a higher back- 
ground level. The observations showed no ob vious additional 



tions from point-like sources. In order to run the source detec- 
tion algorithm of the SAS-software package, we created images 
for the EPIC pn, selecting only single events (PATTERN=0) in 
energy band 1, and single and double events (PATTERN<4) 
for the other bands. For MOS we used single to quadruple 
events (PATTERN<12) in all bands. To avoid differences in 
the background over the EPIC pn detector, we omitted the en- 
ergy range between 7.2 keV and 9.2 keV, where the detector 
background shows st rong spatially variable fluorescence lines 
( IFrevberg et alj|2004[) . We created images, background images 
and exposure maps, and masked them to an acceptable detector 
area. The binning for all images is 2". 

Using the SAS-tasks eboxdetect, version 4.19, and 
emldetect, version 4.60, we searched for point-like sources in 
the field of view (FOV), simultaneously in the 5 energy bands 
and three detectors. First, we searched in each observation sep- 
arately, to correct for inaccuracies in the pointing positions. The 
resulting source lis ts were correlated to catalog ues from USNO 
JMonet et al.ll2003l) . SIMBAeH and Chandra dStrickland et al 



times with solar wind charge exchange (e.g. ISnowden et al.l 120021) . Offsets were determined and applied to each observa 



2004), which would result in times with high background at 



energies below 1 keV, so no further exclusion of exposure time 
was necessary. The exposure times after screening for high 
background (T exp c i ean ) are shown in TableQ] Summing up over 
the final good time intervals, we ended up with 80 ks of expo- 
sure time in total. This means only about 37% of the original 
exposure time could be used for the analysis presented in this 
paper. This number is quite small, compared to typical expo- 
sure time fractions of usable times after screening of 60-709^ 
Next, we screened for bad pixels that were not detected by the 
pipeline. In order to be able to merge images later on, we cal- 
culated sky coordinates (X,Y) for the events in all observations 
with respect to the centre reference position ff2ooo=OO h 47 m 33 s 3, 
<520O0=-25°17'18". For the following analysis we split the data 
set into five energy bands: 0.2-0.5 keV, 0.5-1.0 keV, 1.0- 
2.0 keV, 2.0-4.5 keV and 4.5-12 keV as bands 1 to 5. 

2.1. Point source removal 

In this paper we did not study the population of the point 
sources, but we focused on the diffuse emission in the halo and 
the disc of the galaxy. To do so, we had to remove contribu- 



tion. With the position corrected event files, we again created 
images on which we executed the final point source detection. 
We searched in the merged images from observations 1, 2, and 
4, and separately in the images from observation 3. The rea- 
son for merging only observations 1, 2 and, 4 is that observa- 
tion 3 has a pointing offset (~6') into the northwestern halo and 
therefore we would have different point spread functions on the 
same sky coordinates. 

Additionally, we created a point source catalogue for the 
Chandra observations. Point-like sources in ChandraObslD 
3931 were identified using the Wavelet-Based detection 
Algor ithm (wavdetect in the ciao software, version 3.4, 
Freeman et alj 120021) . in the 0.5-5.0 keV energy band using 
scales of 1", 2", 4", 8", and 16". For ObsID 969 and ObsID 



790 w e adopted the published source list from I Strickland et al 



(120021) . The combined XMM-Newton and Chandra source list 
will be published and further discussed in a forthcoming paper. 
The combined source list was used to remove the point- 
like sources from the data sets. The SAS-task region was used 
to produce elliptical regions that approximate the PSF with an 
analytical model at a given detector position and flux value 
(0.5 times the background flux at this position). Sources that 
were not detected in the XMM-Newton data sets, but are known 



1 see the XMM-Newton EPIC Background Working Group webpage 
http://www.star.le.ac.uk/~amr30/BG/BGTable.html 



http://simbad.u-strasbg.fr/simbad/ 
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Fig. 2. ESO/MPG 2.2m Wide Field Imager (WFI) I-band im- 
age of a possible galaxy cluster in the background of NGC 253. 
The white ellipse marks the source extent of the X-ray source 
as found in the source detection. The bottom quarter of the im- 
age is not covered by the WFI detector. 



from Chandra observations, were excluded with a circular re- 
gion with a diameter of 8". One might argue that these sources 
contribute only little to the overall emission. However, we took 
up a conservative position and also excluded these sources to 
keep any unwanted interference at a minimum. 

One of the sources in the XMM-Newton list in the north- 
western halo at a 2 ooo=OO h 47 m 07 s .7, £ 2 ooo=-25°16'21" (indi- 
cated in Fig. Q]) was flagged as extended by emldetect. In 
a combined I-band image produced from several observations 
with the Wide Field Imager on the 2.2m ESO/MPG telescope 
between July 1999 and August 2000, we also find an enhance- 
ment of sources in this region (see Fig. |2j. These findings in- 
dicate that this source could be a galaxy cluster in the back- 
ground. For the further analysis of the NGC 253 diffuse emis- 
sion the source was removed using a circular region of 1.5' 
diameter. 

2.2. Images 

We used all 4 observations to produce images. The observa- 
tions have different pointing directions and position angles, so 
we obtained images where almost all the CCD gaps are filled. 
The single images, from the energy bands 1 to 5, were corrected 
for the detector background (electronic noise, high energy par- 
ticles) by subtracting the surface brightness of the detector cor- 
ners, that are outside of the field of view. The images were ex- 
posure and vignetting corrected, and adaptively smoothed with 
a Gaussian kernel, with sizes between 10" and 47" (Fig.Q~|i. For 
a detailed description of this method see App.lAl A false-colour 
image was produced by combining the images in the three low- 
est energy bands 1, 2, and 3, as channels red, green, and blue, 
respectively (Fig.0. The image, after the point source removal 
is shown in Fig. [4] 

2.3. Hardness ratio maps and spectra 

As a big advantage, compared to the observations by ROSAT 
and Chandra, the higher count rates in XMM-Newton allowed 




Fig. 3. Adaptively smoothed EPIC pn image of NGC 253. 
The colours correspond to the energy bands (0.2-0.5 keV, red), 
(0.5-1.0 keV, green) , and (1.0-2.0 keV, blue). Overplotted 
in white is the inclination corrected optical D25 ellipse of 
NGC 253. The centre of the galaxy is marked with a black 
cross. Scale and orientation are indicated. 




Fig. 4. Adaptively smoothed EPIC pn image of the diffuse 
emission of NGC 253. Point-like sources have been removed. 
Overplotted in green are the regions that were used for extract- 
ing hardness ratios and spectra. The bright knot in the north- 
western halo (see Sect. 13.2b is indicated by a blue circle. The 
inclination corrected optical D25 ellipse is shown in white. 
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us to extract spectra with reasonable statistics from smaller re- 
gions in the disc and the halo. For the hardness ratios and spec- 
tra we again restricted ourselves to the EPIC pn data. We did 
not use observations 2 and 3 for hardness ratios and spectra, 
because after good time interval screening only little exposure 
was left (cf. TableHJ. 

Energy spectra of several regions (cf. Fig.@]i were extracted 
from the event files, using the SAS-task evselect, version 
3.60, after removal of the point-like sources. To calculate the 
area of these regions, we used the task backscale, version 
1.4.2. We produced background spectra using a region at the 
southwestern border of the FOV, together with observations 
where the filter wheel was closed. A detailed description of this 
procedure, which also handles the binning of the spectra, can be 
found in Appendix IE1 The final, background subtracted source 
spectrum for each region has a significance of at least 3 cr in 
each data bin. The advantage of using this method is a more 
accurate background subtraction which also accounts for vi- 
gnetting and the particle background. As shown in AppendixlBl 
this can have a large effect (of up to ~20%) on the flux values 
and it can eliminate the need of introducing a power law com- 
ponent to account for residuals at high energies in the spectral 
fits. 

Since the emission is mostly confined to energies be- 
tween 0.2 and 2.0 keV, we only calculated the hardness ratios 
HR1 and HR2, where HR1=(B 2 -Bi)/(B 2 +Bi), and HR2=(B 3 - 
B2)/(B3+B2). Bi, B2, and B3 are the count rates in the energy 
bands 1 to 3, i.e. 0.2-0.5 keV, 0.5-1.0 keV, and 1.0-2.0 keV, 
respectively. They were obtained by summing up the back- 
ground subtracted counts in the spectra in the energy bands. 
Observations 1 and 4 were performed for EPIC pn using the 
medium and thin filter, respectively. As these filters have a sig- 
nificantly different throughput, we did not calculate the com- 
bined hardness ratios, but show the results for the two obser- 
vations independently (Fig. [5] and Table |2}. The errors on the 
hardness ratios were propagated from the errors on the counts 
in the spectra. 

In order to fit the spectra with physical models, we created 
the proper response and anxiliary response files for extended 
sources for each spectrum. In XSPEC 11.3.2, we linked the 
model parameters between the two observations and included 
a global renormalisation factor to account for differences be- 
tween the observations to fit the spectra from observations 1 
and 4 simultaneously. To fit the spectra with physical mod- 
els we used the chi-squared statistics in XSPEC. To assess 
the significance of additional model compo nents, the Bayes ian 
Information Criterion (BIC, introduced by ISchwarz 19781 an 
astronomical application can be found, e.g., in iLiddld 120041) 
was applied: a difference of 2 for the BIC is regarded as posi- 
tive evidence, and of 6 or mor e as strong evidence, against the 



mode l with the larger value (iJeffreysl 11961c iMukherjee et al 
1998). We did not use the F-test, since this method is not 



els are the following: (i) the^J is below 1.4, (ii) the BIC differ- 
ence between a given model and a model with more parameters 
is greater than 2, (iii) the errors on all model parameter are con- 
strained, and (iv) the model is physically meaningful. 

3. Results 

To characterise the diffuse emission in the disc and the halo, we 
analysed images in different energy bands, and hardness ratios 
and spectra from several regions. The size of these regions was 
chosen to provide enough counts for a spectral analysis. As for 
the shape, in the disc the regions were chosen in a way that they 
follow the spiral arm structure. In the halo, we chose plane- 
parallel regions above the galactic disc. The projected heights 
of the halo regions are 2 kpc, with exception of the region fur- 
thest to the northwest (region 1), which has a projected height 
of 3 kpc. The regions are overplotted on top the false-colour X- 
ray image in Fig.|4] The hardness ratios in the different regions 
are given in Tableland shown graphically in Fig. [5] 

3.1. Disc diffuse emission 

The disc shows diffuse emission in energies up to ~10 keV, 
where the harder emission is located close to the centre of 
NGC 253. The soft emission (<1 keV) shows the largest extent 
along the major axis. From the nucleus, it reaches -7.0 kpc to 
the northwest and -6.4 kpc to the southeast. 

A prominent feature in the disc is the lack of very soft emis- 
sion northwest of the maj or axis. This is alre ady known from 
ROSAT observations (e.g. IPietsch et al.luOOOl) and can be ex- 
plained by the geometry of the system: The galaxy's disc is 
oriented so that we see the underside of the disc. The emis- 
sion from the northwestern halo behind the disc is therefore 
absorbed by the intervening disc material. 

The spectral properties in different regions of the disc are 
summarised in Table [3] The spectra of all disc regions are 
shown in Fig. [6] A representative example (region 14) of one of 
these disc spectra shows lines from O VII (-0.57 keV), O VIII 
(0.65 keV), FeXVII (0.73-0.83 keV), NelX (-0.91 keV), 
NeX (l.OkeV), MgXI (-1.3 keV), and Si XIII (-1.9 keV). 

To fit the spectra, we tried several different models, which 
all contain an absorption mo del tbabs for the Galactic fore- 
ground JV H of 1 .3 x 10 20 cirT 2 dDickev & Lockmanll990l) . Also 
the a bundances were fix ed to solar values (see also Sec. 14.4b 
from I Wilms et al.l (120001) . A simple one-tem perature thin ther- 
mal plasma model (apec, ISmith et al.l 1200 lb did not result in 



valid to test if an addi tional absorption component is justified 
(IProtassov et al.l 120021) . since the null value of the additional 
parameter is on the boundary of the set of possible parameter 
values, i.e. zero, and negative column de nsity values are not 
supported by the tbabs dWilms et al.l2000l) model. For the BIC, 
these restrictions do not apply. The criteria for our best-fit mod- 



a good fit (i.e. xl ^ 1-4) in any case. Similarly, a power law 
model did not give good fits (see Table [3}. At least three com- 
ponents were necessary for most of the regions: two thin ther- 
mal plasmas plus a power law component. The power law was 
needed to account for the emission above ~ 1 keV and proba- 
bly results from point sources below the point source detection 
limit, or incomplete source removal due to too small extraction 
radii. 

The obtained temperatures are quite uniform throughout 
the disc and vary from 0.1 to 0.3 keV and from 0.3 to 0.9 keV 
for the soft and the hard component, respectively. The intrin- 
sic luminosity (corrected for Galactic absorption) of the diffuse 



M. Bauer et al.: XMM-Newton observations of the diffuse X-ray emission in the starburst galaxy NGC 253 
Table 2. Comparison between hardness ratios, for the individual observations (obs 4: thin filter; obs 1 : medium filter) 



Region 


HR1 


HR2 




Obs 1 


Obs 4 


Obsl 


Obs 4 


Northweste. 


rn halo 








1 


0.16 ±0.06 


-0.04 ± 0.05 


-0.44 ± 0.04 


-0.46 ± 0.04 


2 


0.13 ±0.06 


-0.09 ± 0.05 


-0.41 ± 0.04 


-0.46 ± 0.05 


3 


0.18 ±0.06 


0.10 ±0.04 


-0.50 ± 0.04 


-0.61 ± 0.03 


4 


0.26 ± 0.06 


0.08 ± 0.05 


-0.46 ± 0.04 


-0.51 ±0.04 


5 


0.05 ± 0.05 


-0.06 ± 0.04 


-0.56 ± 0.04 


-0.58 ± 0.03 


6 


0.07 ± 0.04 


0.05 ± 0.03 


-0.52 ± 0.03 


-0.65 ± 0.04 


7 


0.14 ±0.05 


0.12 ±0.04 


-0.54 ± 0.03 


-0.59 ±0.03 


2... 4 


0.19 ±0.03 


0.04 ± 0.03 


-0.46 ± 0.02 


-0.54 ±0.02 


5. ..7 


0.09 ± 0.03 


0.05 ± 0.02 


-0.53 ± 0.02 


-0.62 ± 0.02 


1...7 


0.13 ±0.02 


0.03 ± 0.02 


-0.50 ±0.01 


-0.58 ±0.01 


Disc 


8 


0.38 ± 0.07 


0.39 ± 0.05 


-0.23 ± 0.06 


-0.25 ± 0.05 


9 


0.46 ± 0.06 


0.42 ± 0.05 


-0.31 ±0.05 


-0.30 ±0.04 


10 


0.44 ± 0.08 


0.60 ± 0.07 


-0.18 ±0.06 


-0.14 ±0.05 


11 


0.23 ± 0.05 


0.10 ±0.04 


-0.48 ± 0.04 


-0.54 ±0.03 


12 


0.15 ±0.06 


-0.14 ±0.05 


-0.41 ± 0.04 


-0.60 ± 0.04 


13 


0.43 ± 0.07 


0.47 ± 0.05 


-0.20 ±0.06 


-0.23 ± 0.04 


14 


0.30 ± 0.03 


0.26 ± 0.02 


-0.55 ± 0.03 


-0.57 ± 0.02 


15 


0.01 ± 0.04 


0.04 ± 0.03 


-0.66 ± 0.03 


-0.65 ± 0.03 


16 


0.45 ± 0.04 


0.42 ± 0.03 


-0.17 ±0.04 


-0.28 ± 0.03 


17 


0.31 ±0.03 


0.23 ± 0.02 


-0.45 ± 0.03 


-0.42 ± 0.02 


18 


0.17 ±0.04 


0.12 ±0.03 


-0.55 ± 0.03 


-0.53 ± 0.03 


19 


0.37 ± 0.04 


0.32 ± 0.03 


-0.52 ±0.03 


-0.50 ±0.03 


20 


0.27 ± 0.04 


0.27 ± 0.03 


-0.48 ± 0.03 


-0.39 ±0.03 


21 


0.21 ± 0.06 


0.17 ±0.05 


-0.42 ± 0.04 


-0.51 ±0.03 


Southeastern halo 


22 


-0.02 ± 0.06 


-0.07 ± 0.05 


-0.44 ± 0.04 


-0.58 ± 0.04 


23 


0.02 ± 0.05 


-0.11 ±0.04 


-0.50 ± 0.04 


-0.55 ± 0.04 


24 


0.05 ± 0.08 


-0.12 ±0.07 


-0.25 ± 0.05 


-0.33 ± 0.06 


23+24 


0.03 ± 0.04 


-0.12 ±0.03 


-0.41 ±0.03 


-0.48 ± 0.03 


25 


-0.14 ±0.07 


-0.21 ± 0.06 


-0.26 ± 0.06 


-0.44 ± 0.06 


22. . . 24 


0.01 ± 0.03 


-0.10 ±0.03 


-0.42 ± 0.02 


-0.51 ±0.03 


22. . . 25 


-0.01 ±0.03 


-0.12 ±0.02 


-0.39 ± 0.02 


-0.50 ± 0.02 



emission within the inclination corrected optical D25 ellipse 
is 2.4xl0 39 erg s" 1 (0.2-10.0 keV), or 8.5xl0 38 erg s" 1 (2.0- 
10 keV). Both values were corrected for the area of cut-out 
point-like sources. 

The spectra decrease in hardness from the northwest to the 
southeast parallel to the minor axis, which can easily be seen 
in the hardness ratio maps (HR1, Fig.0. This is not an effect 
caused by different temperatures, but by the increasing strength 
of the soft spectral component towards the southeast (compared 
to the hard component), as the optical depth through the halo 
on the near side of the disc increases. 

For the disc regions that lie in the absorption band (re- 
gion 9) and along the nuclear outflow (regions 16, 17, and 
18) we tested for additional absorption components. In other 
regions we do not expect a measurable absorber in the fore- 
ground with a sufficiently bright source in the background. In 
region 9, we obtain a column density of 0.50*9' "j^xlO 22 cm -2 
for the hard component. The region covering the northwest- 
ern outflow (region 16) closest to the centre shows values of 
0.076+°°^xl0 22 cm- 2 and0.58^^x!0 22 cirr 2 for the soft and 



centre of NGC 253 and the southeastern outflow do not seem 
to require an additional absorption component. 

In a few cases, the two thermal plasma plus a power law 
component model did not give the best fit. The southern re- 
gion on the northeastern end of the disc (region 12) did not re- 
quire a power law component. It does not cover a spiral arm of 
NGC 253, thus the contribution from point-like sources below 
the detection limit may not be significant. The region furthest 
to the south in the disc (region 21) was well fit with one ther- 
mal plasma and a power law. A second thermal plasma was not 
required. 

We found that spectra are harder in regions which cover 
spiral arms. The northeastern regions of the disc (regions 1 1 
& 12) are the best example for this. The region which covers 
the spiral arm (region 11) shows temperatures of 0.18^^ and 
0.58+oJg keV, whereas the region adjacent to the spiral arm (re- 
gion 12) is significantly cooler with temperatures of 0.07 + 0.01 
and 0.25*2™ keV. The latter spectrum is actually more similar 



hard component, respectively. Regions 17 and 18, covering the to the typical halo spectrum (see next section). 
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3.2. Halo diffuse emission 

The halo shows emission only below ~1 keV. Its projected 
maximum extent is ~9.0 kpc to the northwest, and ~6.3 kpc 
to the southeast, perpendicular to the major axis. The gen- 
eral shape resemble s a horn structure. T his was already seen 



with ROSAT (e.g. IPietsch et all l200fi . and Chandra (e.g 



Strickland et alJ I2002J using the term "X-shaped morphol 



ogy"). In the northwestern halo, the EPIC pn images only show 
the eastern horn. In the southeastern halo, both the eastern and 
the western horn are visible in the energies between 0.2 and 
0.5 keV. At higher energies, the western horn is not visible. 

On smaller scales the halo emission seems not to be uni- 
formly distributed. It shows a patchy structure, as was seen be- 
fore in the ROSAT data. One notable feature is a brighter knot 
(see Fig. |4j, which coincides with the nuclear outflow axis in 
the northwestern halo at a height of about 3.5 kpc above the 
disc. It is bright in energies between 0.2 and 1.0 keV and ap- 
pears yellow in the false-colour image (Fig. [3]). We checked 
whether any of the detected structures coincide with chip gaps 
of the detector, and could therefore be artificial, but no correla- 
tion was found. 

The spectral properties in different regions in the halo are 
summarised in Table [4] To fit the spectra, we applied the same 
approach as for the disc. Again, simple models cannot describe 
the spectra. A model with two thin thermal plasmas gave a good 
fit in all regions in the halo. In region 6, and all regions contain- 
ing region 6, adding a power law component results in a better 
fit, possibly pointing at a cosmic ray contribution. 

The southeastern halo is softer than the northwestern halo, 
which results in redder colours in the southeastern halo in the 
EPIC pn false-colour image (Fig.0, and also in lower values in 
HR1. A fit to the spectrum of the whole northwestern halo gave 
temperatures of 0. 10±0.01 and 0.3 1+0.02 keV. The spectrum in 
the southeastern halo is similar, with temperatures of 0.09±0.02 
and 0.29+Q y^ keV. The difference in hardness is because the two 
plasma components contribute different amounts. Compared to 
the normalisation of the hotter plasma, the normalisation of the 
cooler plasma is about 2.1 times stronger in the southeastern 
halo, with respect to the northwestern halo. 

The good statistics of the EPIC pn data allowed a further 
subdivision of the halo into smaller regions. The spectra of all 
halo regions are shown in Figs. [7] and [8] A representative ex- 
ample (region 7) of one of these halo spectra shows prominent 
oxygen lines at 0.57 keV (OVII) and 0.65 keV (OVIII) and 
also an iron line at -0.8 keV (FeXVII). 

The halo is not uniform in its spectral properties on smaller 
scales. The northwestern halo is softer in the east than in the 
west, while the southeastern halo is softer further away from 
the disc (see HR1 maps in Fig. [5}. Additionally, there is a 
marginal indication of a hardening of the emission along the 
direction of the northwestern outflow (regions 3 and 6, signifi- 
cance: I.O-I.60- in HR1 and 0.40-0.80 cr in HR2). 

The total intrinsic luminosity for the diffuse emission, cor- 
rected for the area of the removed point-like sources, in the 
northwestern halo is 8.4xl0 38 erg s -1 (0.2-1.5 keV), compared 
to 2.3xl0 38 erg s -1 in the southeastern halo. To calculate elec- 
tron densities, we assumed a volume for the emitting region. 



We modeled the northwestern halo with a cylinder with a ra- 
dius of 4.3 kpc and a height of 4.2 kpc, plus a cylindrical seg- 
ment with a height of 4.2 kpc, a radius of 4.3 kpc and a width 
in the southeast-northwest direction of 3.0 kpc (to model re- 
gion 1). This gives a volume of 298 kpc 3 or 8.7 x 10 66 cm 3 . 
For the southeastern halo we assumed a cylinder with a radius 
of 3.5 kpc and a height of 2.0 kpc, plus a cylindrical segment 
with a height of 2.0 kpc, a radius of 3.5 kpc and a width in the 
southeast-northwest direction of 3.0 kpc (region 25), resulting 
in a volume of 1 13 kpc 3 or 3.3 x 10 66 cm 3 . To calculate densi- 
ties and the total mass in the emission regions, we corrected the 
volumes for the cut-out point-like sources. Using the emission 
measure of the fit (cf. the documentation of the apec model 
in XSPEC), the resulting densities are 3.2?7 -a5 x 10 -3 cm -3 
and 4.7 jj~ - 5 x 10~ 3 cirT 3 for the northwestern and southeast- 
ern halo, respectively. 77 is t he volume fil l ing fa ctor (77 < 1). 
With solar abundances from IWilms et alJ (120001) . this implies 
total masses of 3.3 77- 05 x 10 7 M Q and 1.8 tT 05 x 10 7 M© for 
the northwestern and southeastern halo, respectively. 
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Fig. 1. Adaptively smoothed EPIC pn images with contours in the lower 4 energy bands: (top-left) 0.2-0.5 keV, (top-right) 0.5- 
1.0 keV, (bottom-left) 1.0-2.0 keV, and (bottom-right) 2.0-4.5 keV. Contours are at (0.35, 0.50, 0.80, 1.6, 2.5, 6.0, 20, 100) x 
10~ 5 ct s" 1 pix -1 . Additionally we show the inclination corrected optical D25 ellipse in black. The black arrow in the lower left 
image indicates the position of a possible background galaxy cluster. 
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Fig. 1. continued, (left): Adaptively smoothed EPIC pn image with contours in the highest different energy band (4.5-12 keV). 
(right): The vignetting corrected exposure map of the merged four observations. The outer contour indicates ks and the exposure 
increases linearly towards the centre by one seventh of the maximum (80 ks) per gray-scale level (0-1 1.4 ks, 1 1.4-22.9 ks, 22.9- 
34.3 ks, 34.3^45.7 ks, 45.7-57.1 ks, 57. 1-68.6 ks, 68.6-80.0 ks). Except for a few pixels, all the detector gaps are covered by at 
least 4.4 ks in the central region. 
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Fig. 5. Hardness ratio maps of observations 1 (left) and 4 (right): The map was binned to the same regions as in Fig. |4] The 
higher the index, the harder the spectrum. The background colour corresponds to an artificially set value. We also show the 
inclination corrected optical D25 ellipse, (top): HR1=(B2-Bi)/(B2+Bi), where Bi and B2 are the count rates in the energy bands 
0.2-0.5 keV and 0.5-1.0 keV, respectively, (bottom): HR2=(B3-B2)/(B3+B2), where B2 and B3 are the count rates in the energy 
bands 0.5-1.0 keV and 1.0-2.0 keV, respectively. 
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Fig. 6. Spectra of different regions in the disc. The red and the black data points and model fits are from observations 1 and 4, 
respectively. The lower panel shows the residuals of the fits. 
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Fig. 7. Spectra of regions in the northwestern halo (as Fig. [6}. 
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Fig. 8. Spectra of different combined regions in the southeastern halo (as Fig. [6}. 



Table 3: Spectral fits to the disk regions. 
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Table 3 - Continued 
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n a q +0.06 
"•-> -0.03 


n. 9+1.0 








1 6 +a4 

'•"-0.3 


1 8+0-3 
1 -8-0.6 


apec + apec 


1.5(176.3/116) 


65.3 




o.i8!8SJ 


3 4 +0 - 5 

J ' -0.9 




64+ 010 

"'" -0.05 


o 4+O.6 
-'■^-0.6 






TO 


apec+apec+po 


0.7(84.5/114) 


-13.4 




o.ins 


9 7+0.8 
^•'-1.0 




59+ 010 
"■•"-on 


9 4+0.8 


1 1+0.6 

'•'-0.3 


1 r 05 




Region 20 
























po 


2.4 (198.2/81) 


81.1 














2 0+ 01 
^•"-o.i 


4 0+° 3 

^•"-0.3 


z 

o 


apec 


3.7(300.1/81) 


115.5 




31+0.01 
"•-' -o.oi 


1.1+g-J x 10 5 












o 


apec+po 


1.1 (83.0/79) 


17.7 




0.30!°;* 


4.0+|-° 








1 7+ 03 

'•'-0.2 


9 c+0.6 
z — '-0.5 


IO 


apec + apec 


2.6 (205.9/79) 


93.1 




0-i8!SS 


4 0+ 16 
t.u _ 05 




64+ 8512 

"*" -0.08 


q 7+0.5 
"'■'-2.1 








apec+apec+po 


1.0(73.7/77) 


16.6 




0.21+™ 5 


3 +LI 

-'■"-1.8 




57 +0 16 

"■-"-0.47 


1 7 +as 

1, '-1.6 


1 5 +a4 

1 -'-03 


9 9+0.7 
^ -0.7 





Continued on next page. 



Table 3 - Continued 



Region/model 



xiix'M 



BIC 



(10 22 cm" 2 ) 



kT 

(keV) 



norm 

(io- 5 ) 



(10 22 cm" 2 ) 



kT 
(kev) 



norm 

(io- 5 ) 



(io- 5 ) 



Region 21 
po 
apec 
apec+po 

apec + apec 
apec+apec+po 



1.2(48.8/41) 


13.0 


2.1 (87.7/41) 


38.2 


0.8 (32.3/39) 


2.7 


1.1(43.2/39) 


15.2 


0.8 (30.2/37) 


7.4 



26 +003 
20+ 007 

"■■'"-0.07 

0.16^ 
0.18!^ 



l.l + °- 2 xl0 5 
1 7 +0 - 7 
3.3 + '-° 



1.6 



-0.8 
1.1 
1.6 



0-63^ 



1 6 +0 - 4 
4S+ - 65 



1 7 +02 



1 4+0-6 



1 5 +07 



T 8+0-5 



9 i +0.8 
' i,1 -0.7 



1 7 +1 2 



n.c: not constrained. 

NOTE: All errors are 90% confidence for a number of interesting parameters equal to the number of free parameters in the model (2, 4, 6, 7, 8, or 9 free parameters depending on the 
model). Thin thermal plasma models are apec, absorption models are tbabs. All models have a fixed foreground absorption of 1.3 x 10 20 cnr 2 . The preferred model is printed in bold font. 



ffi 



X 

z 



o 
cr 



5 



X 



p 

3- 



z 

o 
o 
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Table 4: Spectral fits to the northwestern halo regions. 



Region/model 


xi(x 2 /y) 


BIC 


kT 


norm 


kT 


norm 


r 


norm 








(keV) 


(io- 5 ) 


(kev) 


(IO" 5 ) 




(io- 5 ) 


Northwestern halo 


Region 1 


















po 


1.1(25.9/23) 


7.3 










2 0+°- 3 

z '"-0.3 


6 1 +M 


apec 


3.9 (88.6/23) 


38.1 


25 + 002 

" -0.02 


i.r°- 3 x io 5 










apec+po 


1.1 (22.1/21) 


9.8 


19+ 021 

"• 1? -0.19 


1 3+ 29 

1 -°-1.3 






i q+0.6 
l - y -0.6 


5 6+ 18 

J *"-2.0 


apec+apec 


1.0(21.1/21) 


8.6 


10+° ° 2 

" x -0.02 


19+ 8.5 

17 -8.3 


35+ 006 

-0.05 


5 6+ L6 

•'•"-1.5 






apec + apec+po 


0.9(16.3/19) 


8.6 


09+° ° 2 

"•" 7 -0.07 


20+ 21 
z.o_ 20 


29+ 020 

-0.29 


2 6+ 31 

^•"-2.6 


55+ 175 


c q+4.0 
J.J_ 36 


Region 2 


















po 


1.1(21.7/20) 


5.9 










2 0+°- 4 

^■"-0.4 


q 1+0.9 

- > - 1 -0.8 


apec 


2.6(51.4/20) 


24.9 


24+ 002 

"•^ -0.02 


i.i + °ix io 5 










apec+po 


0.7 (12.5/18) 


-0.1 


n O4+0.08 

"-^-0.06 


™ + -\i 






9 r+0.6 
^ -0.9 


1.5+- 


apec+apec 


0.8 (14.6/18) 


3.3 


10+° ° 3 

"- 1 -0.02 


8-3!« 


O Q9+0.06 
"•-' -0.05 


3 0+ 1 ° 

-'•"-0.9 






apec + apec+po 


0.7(11.8/16) 


4.8 


16+ nxJ 
u.iu _ 0)6 


1 4+2.5 


031+ n.c.' 

"■-"-031 


1 7+ 2 ' 1 


9 Q+0.9 


0.84+- 4 


Region 3 


















po 


1.6(67.0/41) 


26.6 










1 7+ 03 

1,1 -03 


4 5+0.8 
^-0.7 


apec 


2.4 (97.9/41) 


42.9 


29+° ° 2 

"•^ -0.02 


l.l+°- 2 x 10 5 










apec+po 


1.1 (41.2/39) 


13.2 


O ao+0.06 
"■-' -0.04 


O Q+1.4 

3.8_ u 






3 r°' 3 

-'■ 1 -0.6 


0.87!° 54 


apec+apec 


0.9 (35.3/39) 


6.6 


09+° 01 

"•" -0.02 


12+ij 


34+ 005 

"•-' -0.03 


4 6 + ° 9 

t-"-0.9 






apec + apec+po 


0.9 (34.0/37) 


12.5 


09+° ° 2 

"•" 7 -0.02 


121 


33+ 007 

U,J -0.05 


4 2+13 


-0.1 1+ 277 , 
-n.c. 1 


1 o+ 2 -' 

x -"-i.o 


Region 4 


















po 


1.7(45.1/27) 


19.5 










1.8+°; 3 


q 9+0.7 
-'■ z --0.6 


apec 


1.8(47.8/27) 


21.2 


26+° ° 2 

"•^"-0.02 


1.2 + °f X 10 5 










apec+po 


1.0(24.4/25) 


8.5 


29+° ° 5 

"""-0.05 


3 r 13 

J - -1.3 






9 Q+0.6 
^■'-l.l 


71+0.59 
"• ' 1 -0.68 


apec+apec 


0.7 (16.5/25) 


-2.9 


11 + 004 

"- 1 -0.03 


5 T 3 - 9 

J -2.9 


33+ 007 


■2 a +0.9 
-3-0.9 






apec+apec+po 


0.7(16.4/23) 


3.7 


11 + 008 
"• xl -0.04 


5 6+ 40 


34+ 009 

"•-' -0.05 


*2 9+1.4 
-'•^-1.2 


-0.96--; 


0.0+ 28 


Region 5 


















po 


1.5 (63.3/43) 


23.0 










9 q+0.2 
z -0.2 


2.8!" 


apec 


2.5(109.3/43) 


47.5 


23+ 001 

"• z '-'-0.03 


i.r°- 2 x io 5 










apec+po 


0.8 (32.9/41) 


1.1 


21+ 004 

u ' z ' 1 -0.04 


2 5+ 12 






2 6+°- 4 

^•"-0.5 


1 5+ 06 


apec+apec 


0.8(31.1/41) 


-1.4 


0.10+™ 2 


12+ 4 


30+ 004 

"■•'"-0.04 


T 4+O.6 
J - -0.9 






apec+apec+po 


0.7 (28.8/39) 


2.8 


0.081°$ 


13+10 

10 -n 


23+ 009 

"•"-0.06 


9 7+1.4 
''■'-1.4 


1.6+" . 
-n.c. 1 


1 6+ 15 

x -"-1.6 


Region 6 


















po 


1.5 (141.9/92) 


47.8 










9 9+0.1 

^•^-o.i 


4 6+ 05 

t.u _ 05 


apec 


2.5 (233.0/92) 


94.4 


25 + 001 


i.o!°-; x io 5 










apec+po 


0.8 (74.9/90) 


-3.2 


25 + 003 
"•■"-0.03 


4 2+14 






9 4 +0.3 


2 6+ 06 

^•"-0.7 


apec+apec 


0.9 (82.5/90) 


5.9 


10+ 001 

" -0.01 


15+ 3 

1J -4 


32+ 003 

"•-' -0.02 


6 0+° 9 

"•"-0.8 






apec+apec+po 


0.7 (65.7/88) 


-6.4 


09+° ° 2 

"•" 7 -0.02 


13+J 


a 9Q+0.04 
"■■"-0.04 


4 5+1.4 


1 6+°- 8 

-2.9 


2 1+ 15 


Region 7 


















po 


1.2 (70.9/61) 


15.7 










1 S+ - 2 


5 l+°- 8 


apec 


1.7(103.9/61) 


39.8 


25 + 001 


i.i!°- 2 x io 5 










apec+po 


0.7(41.7/59) 


-9.4 


n 99+0.04 
u - zz -0.04 


9 q+18 
^• 7 -1.5 






1 7+ 05 

x ''-0.6 


3 5+ 11 

-'•-'-1.5 


apec+apec 


0.9 (52.0/59) 


4.5 


11 + 002 
"• xl -0.03 


8-2!^ 


Q3J+0.03 


5 0+° - 9 
J -"-i.i 






apec+apec+po 


0.6 (33.5/57) 


-14.9 


06+° 01 

"•""-0.02 


25+ 53 


22+ - 05 

"""-0.05 


T T+2.6 

•'•■'-1.6 


0.70+ 110 , 
-n.c. 1 


4.2+^ 


Regions 2. . . 4 


















po 


1.5 (147.2/101) 


46.0 










2 0+°- 2 

^■"-0.2 


9.8+|; 2 


apec 


2.2 (219.7/101) 


87.3 


0.26!°,;°,! 


l.l + °-J x IO 5 










apec+po 


1.0(101.6/99) 


17.1 


0.28!°« 


7 5+ 2 - 9 






2 5+°- 3 


4 5+ 14 


apec+apec 


1.0(97.8/99) 


13.2 


10+ 001 

" -1 -0.01 


27+9 


33+ 003 

"• J -0.03 


H-1 






apec+apec+po 


0.9 (90.5/97) 


14.5 


09+° ° 2 

"•" 7 -0.02 


25I} 1 


031+ o.o5 

"""-0.04 


9-0!^ 


1.2+ 12 . 
-n.c. 1 


3 5+ 4 ' 1 

-'-'-3.4 


Regions 5. . .7 


















po 


2.0(341.4/170) 


128.2 










9 4+0.1 


10!| 


apec 


2.3 (387.5/170) 


150.0 


25 + 001 

"•^ -0.01 


l.l + °-J x 10 5 










apec+po 


0.9(151.8/168) 


-0.9 


25 + 002 

"• -0.02 


12+j 






9 7+O.2 
^•'-0.2 


4.8!| : ° 


apec+apec 


0.9(157.8/168) 


5.8 


10+ 001 
"■ lu -o.oi 


37+* 


32+ - 02 

"•-' -0.02 


15+J 






apec+apec+po 


0.8(137.8/166) 


-7.2 


10+° ° 2 

"• lu -0.02 


28+12 


29+ 003 

"•^'-0.03 


12+ 3 


1 9+°" 7 


3.2+ 2 ; 2 



Continued on next page. 



M. Bauer et al.: XMM-Newton observations of the diffuse X-ray emission in the starburst galaxy NGC 253 



17 



Table 4 - Continued 



Region/model 


xi(x 2 M 


BIC 


kT 


norm 


kT 


norm 


r 


norm 








(keV) 


(io- 5 ) 


(kev) 


(IO" 5 ) 




(io- 5 ) 


Regions 1 ... 7 


















po 


2.2 (576.0/263) 


216.9 










2-8_2! 


22+1 


apec 


2.7 (720.2/263) 


276.1 


25+°°' 

"■"-0.01 


1.1+21 x 1Q5 










apec+po 


0.8 (201.6/261) 


-50.1 


26+° 01 
"■^--o.oi 


29+ 4 






9 O+0.2 
Z -°-0.2 


IO! 2 


apec + apec 


0.8 (209.4/261) 


-40.1 


10+°°' 
"- 1 -0.01 


98!|? 


33+ 002 
"• J -o.oi 


35+ 3 






apec+apec+po 


0.7 (170.2/259) 


-83.9 


10+°°' 
"■ lu -o.oi 


64+ 14 


031 +0.02 
J -0.02 


29+ 5 


9 4+0.4 
^ M).6 


6.6-l : ; 


Southeastern halo 


Region 22 


















po 


1.3(64.1/50) 


18.8 










2 5+°- 3 


i a +o.4 
i-8_o. 4 


apec 


3.1 (154.0/50) 


64.4 


12+° 13 
"• lz --o.oi 


1.2+2 3 xl0 5 










apec+po 


0.9 (44.0/48) 


7.1 


16+° ° 6 

"- 1 -0.05 


1-8-11 






9 4+0-6 


1<1 


apec+apec 


0.9 (45.3/48) 


8.6 


10+° ° 2 

"- 1 -0.02 


8.9!" 


n Q9+O.O6 
"•-' z --0.05 


1-9-21 






apec+apec+po 


0.8 (35.5/46) 


3.9 


07+° ° 4 

"•"'-0.02 


15+ 16 


1 8+0.09 


i.8*l:! 


0.68-- , 


1.8+li 


Region 23 


















po 


1.7(108.7/64) 


41.3 










9 8+0-2 


1.8+2 3 


apec 


2.6 (168.2/64) 


70.1 


23+ 001 
"■■"-o.oi 


i.i+2 2 x io 5 










apec+po 


0.9 (52.8/62) 


2.0 


24+° 03 

"•^ -0.03 


9 7+1.0 
^•'-0.9 






t -3+0.4 
-'-0.4 


o.80!2 : I 


apec+apec 


0.9 (54.0/62) 


3.5 


09+° ° 2 

"■" -0.01 


u:i 


30+ 004 

"■ -0.03 


3 l+°- 8 

J - -0.7 






apec+apec+po 


0.8 (49.2/60) 


5.7 


07+° ° 2 

"•"'-0.02 


17+" 

1 '-10 


25+°"° 5 

"•"-0.04 


3 0+'-° 


2 5+°" 9 


n £7+0.72 
0.67_ 067 


Region 23+24 


















po 


1.5(114.0/77) 


37.7 










2 5+ 02 

•^ -0.2 


t q+0.7 
-^-o.e 


apec 


4.0(309.1/77) 


116.5 


12+° " 
"• lz --o.oi 


i.i+2 2 x io 5 










apec+po 


1.3(98.6/75) 


35.0 


22+° 07 

"■"•-0.08 


91+ 4 - 72 

"■ ; ' 1 -0.50 






2 5+°- 3 


3 5+ 03 

J - -1.9 


apec+apec 


1.1 (82.5/75) 


20.9 


07+° 01 
"■"'-o.oi 


47+21 

H '-15 


26+ 003 

""'-0.02 


6 2+ L0 

"■^-1.0 






apec+apec+po 


0.8 (62.0/73) 


7.1 


06+° 01 
"•""-o.oi 


88+} 16 


21+ 005 

"• zl -0.02 


4 0+ 2 3 

H -"-1.9 


1.1 + 13 . 
-n.c. 1 


3 r°- 8 

->-2.0 


Region 25 


















po 


1.6 (35.8/22) 


16.0 










3 l+°- 5 

-'■ 1 -0.4 


0.83! 027 


apec 


1.2(27.1/22) 


9.3 


o.i2*22i 


i.o+2 3 x io 5 










apec+po 


0.5(10.4/20) 


-7.4 


0.18+^ 


1-9^ 






o 8+0.7 


27+ 018 

"•"-0.26 


apec+apec 


0.4 (8.5/20) 


-12.2 


0.09+22 3 


8.5+ 33 


r, 94+O.O8 
"•^ -0.06 


1 3+ L2 

-0.8 






apec+apec+po 


0.5 (8.4/18) 


-6.1 


0.08tJJ 


8.5!™ 


20+ 007 

"•^"-0.05 


1 7+ 2J 
1 -'-l.l 


2.9+ 18 . 
-n.c. 1 


14+ 0J4 

"• 1 ^-0.14 


Regions 22. . . 24 


















po 


1.9(158.7/83) 


62.0 










9 7+O.2 
"'-0.2 


4 1+0.6 


apec 


3.9(319.6/83) 


121.5 


12+° 16 
"■ lz --o.oi 


1.2+2^ x 1q5 










apec+po 


1.2(93.9/81) 


26.2 


23+° 03 

"■"-0.03 


5 7+ 20 
j./_ 18 






3 r°- 3 


i q+0.6 
x - ? -0.6 


apec+apec 


1.1(91.1/81) 


23.7 


09+° ° 2 

"•" -0.02 


29+ 7 


29+°-° 3 

"•"-0.03 


6 7+ 17 

"•'-1.4 






apec+apec+po 


1.1 (83.5/79) 


25.1 


07+° 01 
"•"'-o.oi 


49+29 

*'-22 


24+ 003 

"•^-0.04 


6.6! 2 ;° 


1.7+ 12 . 
-n.c. 1 


1 7+ L8 

1 -'-1.5 


Regions 22. . . 25 


















po 


1.7(130.9/78) 


48.2 










9 7+O.2 
"'-0.2 


5 4+°- 9 

J - -0.9 


apec 


2.7 (213.8/78) 


87.4 


22+° 01 

"■"•-0.03 


1.2+2 2 x 10 5 










apec+po 


1.1 (81.2/76) 


18.7 


21+° 04 

"■^ -0.03 


6-2^ 






3 0+°- 4 

J -"-0.4 


9 7+O.8 

^■'-o.s 


apec+apec 


1.1 (82.4/76) 


19.9 


09+° ° 2 

"■"-'-0.02 


35+' 


29+ 003 

"•"-0.04 


7 5+ 20 

'• J -1.8 






apec+apec+po 


1.0(70.7/74) 


16.4 


07+° 01 
"•"'-o.oi 


661H 


q 91+0.05 

"• zl -0.04 


7 0+ 2 7 

'•"-2.9 


1.6+ 12 , 
-n.c. 1 


3-2!^ 



n.c: not constrained. 

NOTE: All errors are 90% confidence for a number of interesting parameters equal to the number 
4, or 6 free parameters depending on the model). Thin thermal plasma models are apec, absorption 
a fixed foreground absorption of 1.3 x 10 2 ° cirr 2 . The preferred model is printed in bold font. 



of free parameters in the model (2, 
models are tbabs. All models have 
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Table 5. Flux and luminosity values in the extraction regions. 



Region 



Fx 
(erg s~' cnr 2 ) 



(ergs" 1 ) 



Northwestern halo 










1 


1.7^ x 


10" 


3 (0.2-1.0 keV) 


1.4 x 10 38 


2 


8.4 +2 - 7 x 


10" 


4 (0.2-1.0 keV) 


6.7 x 10 37 


3 


1 1 + 03 X 
1 - 1 -0.4 X 


10" 


3 (0.2-1.0 keV) 


8.8 x 10 37 


4 


8.7 + f 3 x 


10" 


4 (0.2-1.0 keV) 


6.9 x 10 37 


5 


1.0* 3 x 


10" 


3 (0.2-1.0 keV) 


8.0 x 10 37 


6 


1 5 +0 - 4 x 
1.3 _„ 2 X 


10" 


3 (0.2-1.0 keV) 


1.2 x 10 38 


7 


1.2!- x 


10" 


3 (0.2-1.0 keV) 


9.6 x 10 37 


2... 4 


2.8$ x 


10" 


3 (0.2-1.0 keV) 


2.2 x 10 38 


5.. .7 


4.4*J x 


10" 


3 (0.2-2.0 keV) 


3.5 x 10 38 


1...7 


1.0*1 x 


10" 


2 (0.2-2.0 keV) 


7.9 x 10 38 


Disc 




9 


8 - 2 -Z8 X 


10" 


4 (0.2-4.0 keV) 


6.5 x 10 37 


11 


1.3*«x 


10" 


3 (0.2-2.0 keV) 


1.0 x 10 38 


12 


6.9!- x 


10" 


4 (0.2-1.0 keV) 


5.5 x 10 37 


13 


7.2^x 


10" 


4 (0.2-3.0 keV) 


5.7 x 10 37 


14 


2.5 02 x 


10" 


3 (0.2^1.0 keV) 


2.0 x 10 38 


15 


1.5!- x 


10" 


3 (0.2-2.0 keV) 


1.2 x 10 38 


16 (NW outflow) 


1.9*| x 


10" 


3 (0.2-5.0 keV) 


1.5 x 10 38 


17 (Centre) 


2.6!^ x 


10" 


3 (0.2-5.0 keV) 


2.1 x 10 38 


18 (SE outflow) 


1.5!- x 


10" 


3 (0.2-3.0 keV) 


1.2 x 10 38 


19 


1 5+— x 
1.3 _ 02 X 


10" 


3 (0.2^1.0 keV) 


1.2 x 10 38 


20 


2-0!°i x 


10" 


3 (0.2-5.0 keV) 


1.6 x 10 38 


21 


6.0!li x 


10" 


4 (0.2-1.0 keV) 


4.8 x 10 37 


Southeastern halo 




22 


7.7«* x 


io- 


4 (0.2-1.0 keV) 


6.1 x 10 37 


23 


1.0!" x 


10" 


3 (0.2-2.0 keV) 


8.0 x 10 37 


23+24 


1.7!- x 


10" 


3 (0.2-1.0 keV) 


1.4 x 10 38 


25 


5 5+ 2 - 2 x 
j.j _ 30 A 


10" 


4 (0.2-1.0 keV) 


4.4 x 10 37 


22. . . 24 


2.3_ ; 7 x 


10" 


3 (0.2-2.0 keV) 


1.8 x 10 38 


22. . . 25 


2.7*| x 


10" 


3 (0.2-2.0 keV) 


2.2 x 10 38 



4. Discussion 

4.1. The extent of the diffuse emission of NGC 253 

Extended emission from the soft no rthwestern halo was first re- 
ported from Einstein observations dFabbianoll 1988k . Later, ob- 
servations with ROSAT also dis covered the southeastern halo in 
X-rays (e.g. lPietsch et al.l2000k . The ROSAT images in the soft 
band trace the emission in the outer halo to projected distances 
of up to 9 kpc, both in the northwest and the southeast direction. 
With XMM-Newton, the emission is detected out to 9.0 kpc to 
the northwest and 6.3 kpc to the southeast. This difference in 
the southeastern halo can be explained by the high ROSAT sen- 
sitivity extending down to 0. 1 keV. The useful XMM-Newton 
EPIC pn range is limited to 0.2 keV. This makes a big differ- 
ence as there are many strong lines from O IV, Ne VIII, Mg IX, 
MgX, Si IX, and SiX in the energy band between 0.1 and 
0.2 keV For a thermal plasma at a temperature of ~0.1 keV, 
these lines are even stronger than the O VII and O VIII lines, 
and about 60% of the total flux in the energy band from 0. 1 to 
2.0 keV originates from lines below 0.2 keV The southeastern 
halo shows softer emission than the northwestern halo, there- 
fore, the effect is strongest in the southeastern halo. 



Also in the disc the extent of the emission is different. The 
ROSAT images trace the soft emission ~6.8 kpc towards the 
northeast and ~5.3 kpc towards the southwest. With XMM- 
Newton, the disc emission has an extent of ~7.2 kpc and 
~6.3 kpc to the northeast and southwest, respectively. The 
disc spectra are harder than the halo spectra, and therefore the 
higher XMM-Newton sensitivity at energies >0.4 keV comes 
into play. 

4.2. Is the diffuse emission in the disc really from hot 
interstellar gas? 

In this section we investigate the possibility that the observed 
diffuse emission originates in fact not from a truly diffuse hot 
interstellar gas component, but from a large population of weak 
stellar-type X-ray sources that shows the same characteristics. 
This effect was fi rst discovered in the M ilky Way's ridge 
X-ray emission (e.g. iRevnivtsev et al.l 120061) . where evidence 
was found that the bulk of the Galactic ridge X-ray emission 
is composed of weak X-ray sources, mostly cataclysmic vari- 
ables and coronally active binary systems, with a luminosity 
of most of these sources of less than 10 3 erg s" 1 . Also in other 



galaxies. lRevnivtsev et al.l (120071) found that the apparently dif- 
fuse emission is consistent with the emission from an old stellar 
population like in the Milky Way. 



Following the method by IRevnivtsev et al.l (120071) . we used 
K-band observations, to infer the emissivity of the diffuse X- 
ray component per unit stellar mass. We derived the near- 
infrared luminosity and stellar m ass of NGC 253, using the 
total K-band magnitude of 3.772 djarrett et al.ll2.Q03b . the dis- 
tance modulus of 27.06, co rrected for interstellar extinction of 
0.007 dSchlegel et al.ll 1998k. and the colour-de pendent K-band 
mass-to-light ratio from lBell & de Jona d200lk , log(M„/L K ) - 
-0.69 2 + 0.652 x(B-V), with (B-V) = -0.16JCo mer6net al 
2003k . This yielded a total K-band luminosity Lk = 1-6 x 
10 n L Q and a total stellar mass M« = 2.6 x 10 10 M Q . With a X- 
ray luminosity of NGC 253 of 2.0xl0 39 erg s" 1 (0.5-10 keV), 
the emissivity of the diffuse X-ray component per unit stellar 
mass then resulted in 



^0.5-lOkeV 

M, 



7.6!^(±2.3)xl0 28 ergs" 



MZ 



The errors are statistical errors on the measured X-ray flux. 
Additionally, we assumed an uncertainty of ~30% (given in 
parenthe ses), which migh t be associated with the Lk to M t con- 
version dBell et alj|2003l) . The emissivity of NGC 253 should 
only be considered a lower limit. We cut out a quite large re- 
gion in the centre of NGC 253 and corrected for this by filling 
the hole with the average flux of the disc. Therefore, the ob- 
tained X-ray luminosity as well as the emissivity are probably 
underestimated. 

From the luminosity an d other properties of th e Galactic 
ridge X-ray emission (e.g. IRevnivtsev et al.l 120061) and from 
direct measurements of t he luminosity functio n of sources in 
the solar neighbourhood dSazonov et al.l 120061) . the combined 
0.5-10 keV emissivity of cataclysmic variables and coronally 
active stars has been estimated as Lx/M* ~ 1.2 ± 0.3 x 
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Fig. 9. Comparison between the K-band (left, 2MASS. iJarrett et al.ll2003l) and X-ray (right) morphology. The white contours 
overplotted on the X-ray image represent the K-band brightness levels. Both images are on a linear colour scale, and on the same 
spatial scale as indicated in the K-band image. 



10 erg s Mq , The value derived for NGC 253 is larger 
than the value for the Milky Way, indicating the presence of 
a hot gaseous component. 

An even stronger argument is the following: if the diffuse 
X-ray emission is produced by an old stellar population, then 
their morphologies should be similar. A compa rison of the 



X-ray emission with the 2MASS K-band image (IJarrett et al 



20031) is shown in Fig. [9] We found that the X-ray morphology 
does not match the K-band morphology, therefore the diffuse 
emission is indeed not simply due to an old stellar population, 
but has to have a truly diffuse component. 

4.3. Spectral fits and variations in the halo 

As it was mentioned already in earlier publications, there is an 
ambiguity in the spectral fits between a pure multi-temperature 
thermal plasma model and a combinati on of therma l plas- 
mas plus a power la w component (e.g. iDahlem et alj 120001: 
Strickland et al.l2002l) . This ambiguity in the halo emission still 
exists with the XMM-Newton data. Fits to the halo spectra with 
a thermal plasma plus a power law model (see Table|4]i resulted 
in similarly, as for a multi-temperature thermal plasma model. 
A power law component from point sources could be excluded, 
since we were careful to remove any point source contribution. 
Another source for non-thermal emission could be synchrotron 
emission from cosmic ray electrons that are advected with the 
superwind or are accelerated locally in internal wind shocks. 
This is possibly what we observe in the spectrum of region 6 
(and spectra containing region 6). This region lies along the 
northwestern outflow axis. Here the spectral fit improves sig- 
nificantly by adding a power law component to the two thermal 
plasmas. A comparison of the X-ray emission to the 330 MHz 
and 1.4 GHz radio emission (ICarilli et al.lll992l) showed that 



the radio emission is more extended, and does not show the 
horn structure that we see in X-rays. Because of this inequality, 
we prefer the multi-temperature thermal plasma over the single 
temperature model plus power law model for most of the X-ray 
halo emission at the moment. 



The appearance of a "multi-temperature" halo is a natu- 
ral consequence of the plasma driven out of collisional ion- 
isation equilibrium (CIE) towards a non-equilibrium ionisa- 
tion (NEI) state in a galactic outflow. It has to be borne in 
mind that the processes of ionisation, recombination and adi- 
abatic cooling occur all on different time scales. Whereas the 
former two are controlled by microphysics the latter is due 
to t he hydrodynamics of the flow. As h as been described by 
e.g. lBreitschwerdt & Schmutzlerl (119941) adiabatic cooling can 
become the shortest in an accelerating outflow, like the one 
occurring in a starburst driven wind, leading to delayed re- 
combination of the plasma. This means essentially the pres- 
ence of highly ionised species in a plasma flow with low 
kinetic electron temperature, similar to the ions in the so- 
lar wind. Releasing the energy stored in these highly ionised 
species results in diffuse soft X-ray emission and some spec- 
tral hardening. Observing this with the limited spectral res- 
olution of Chandra or XMM-Newton CCDs, reveals a multi- 
temperature halo, although physically the plasma does by no 
means have "different temperatures", but is characterised by a 
non-equilibrium distribution of ions, which by itself is the re- 
sult of the thermodynamic path the system has taken and re- 
flects directly the hydrodynamic flow characteristics. A cur- 
rently ongoing analysis with non -equilibrium models (e.g. 
Breitschwerdt & Schmutzlerl 1 1999b might also be able to ex- 
plain the observations and is the subject of a forthcoming pa- 
per. 
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The northwestern halo shows significant hardness varia 



tions i n HR1, as opposed to the findings by IStrickland et al 



(2002). We checked if this can be caused by a different energy 
band selection, but the result is indep endent whether we use 
the bands from lStricklandet all d2002h (0.3-0.6 keV and 0.6- 
1 .0 keV) or our own. These hardness variations might also be a 
sign of NEI X-ray emission. 

4.4. Temperatures, abundances and column densities 

The X-ray emission from NGC 253 has been observed before 
with several other X-ray observatories. Especially the early ob- 
servatories did not allow to separate the point sources from 
the diffuse emission since the point spread function was quite 
large. Hence, only a combined fit of the emission from the halo, 
the disc, and the nuclear region was possible. Temperatures 
of multi-temperature models ranged between 0.1 and 0.3 keV 
for the low, and bet ween 0.6 and 0.7 keV for the high tem- 



perature componen t JDahlem et alj Il998c IWeaver et al.l uOOOt 
Dahlem et al.ll2000l) . Reported abundances were mostly highly 



subsolar and therefore unphysical for a supposedly metal en- 
riched starburst galaxy plasma. 

Only the X-ray observatories like XMM-Newton, Chandra, 
and to some degree ROSAT allow us to separate the halo from 
the disc emission and to remove contribut ion by point sources 
via a spatial selection. From ROSAT data. lPietsch et al.l (120001) 
inferred a foreground absorbed two-temperature thermal model 
with temperatures of 0.13 and 0.62 keV for the northwestern 
halo emission. No highly subsolar abundances were required. 
The disc emission could be explained by a 0.7 keV thermal 
plasma and an additional thermal plasma (kT=0.2 keV) in front 
of the disc coming from a coronal component. 



The first XMM-Newton results by iPietsch et al.l (120011) re- 
quired a two-component model for the disc emission with tem- 
peratures of 0.13 and 0.4 keV plus residual harder emission, 
possibly from unresolved point sources. For the nuclear region 
three temperatures were needed (0.6, 0.9, and 6 keV). Both 
models used solar abundances. No analysis of the halo emis- 
sion was presented in their paper. 

The best spatial resolution is provided by the Chandra 
observatory. Results on the diffuse dis c and halo emission 
were first published by IStrickland et al.l (12002b . For the halo 
emission they needed a multi-temperature model (apec) with 
at least two temperatures of 0.24 and 0.71 keV (the latter 
with quite large errors) and with a foreground absorption of 
5.3 X 10 20 cirT 2 . A power law (F = 3.3) plus a thermal 
plasma (kT=0.24 keV) gave a similarly good fit. A combina- 
tion of other thermal models (mekal) or non-equilibrium mod- 
els (vnei) did not result in better fits. The diffuse emission from 
the disc was fitted with the same models, however, the temper- 
atures were lower than in the halo, with 0. 17 and 0.56 keV, re- 
spectively. The foreground absorption yields 4.7x 10 20 cm~ 2 . In 
all cases unphysically sub-solar abundances had to be assumed. 

The temperature values for the halo emission, as found by 
our analysis, are lower than the ones from previous observa- 
tions. Our soft component is about 0. 10 keV, which is still com- 
patible with the ROSAT results. However, the hard component 



is only ~0.31 keV for the northwestern halo and ~0.29 keV 
for the southeastern halo. A higher temperature was not neces- 
sary in any of our fits. A possible explanation could be the way 
the spectra were background subtracted. We used a sophisti- 
cated method (see App.lBl that uses the local background at the 
border of the field of view, where no emission from NGC 253 
is ex pected, while other aut hors used e.g. blank-sky observa- 
tions (IStrickland et al.ll2002l) . Using a background from differ- 
ent times and different fields on the sky can lead to systematic 
effects in the background substraction. A background region 
with a higher contribution of the local bubble could, for exam- 
ple, lead to an over-correction, especially at very soft energies 
(<0.5 keV). 

In the disc we found temperatures between 0. 1 and 0.3 keV 
and between 0.3 and 0.9 keV, for the soft and the hard compo- 
nent, respectively. This is consistent with earlier results. 

We also tried to constrain the metal abundances in our fits. 
However the errors on the obtained values are so large, that we 
are not able to give well constrained abundances (northwest- 
ern halo: Z = 0.3^2 Z Q , southeastern halo: Z = 0.4^3 Z , 
disc: Z = 1.0^2'? Z ). Since we do not expect highly subso- 
lar abundances in an environment which is enriched with met- 
als from the starburst via the superwind and galactic fountains, 
we fixed the abundances in our analysis to solar. This is en- 
tirely consistent with the above values. A reason for the low 
abundances, found with different instruments, could be that due 
to a lower spatial resolution and or/and sensitivity more point 
sources contribute to the final spectrum, increasing the contin- 
uum flux. The ratio of line emission to continuum flux is there- 
fore decreased, which mimics the spectral shape of a plasma 
with low metal abundances. A similar effect can be achieved 
when a NEI spectrum is fitted with CIE models. Also, a too 
simplistic model could be the reason, combining regions with 
different temperatures. 

In the disc, two of the regions required an extra absorption 
component in the spectral model. The additional column densi- 
ties range between 0.76xl0 21 cm -2 and 5.8xl0 21 crrT 2 . Direct 
radio measurements of the H I column density showed lower 
values than we derived from the X- ray data. An interpolation of 
the H I maps bv lPuche et all dl99ll) and lKoribalski et al.l d 19951) 
resulted in (~2.4 and 3.4)xl0 21 cirT 2 for regions 9 and 16, re- 
spectively. However, the H I value for the region including the 
nucleus of NGC 253 (region 16) is affected by HI absorption, 
so the resulting column density can only be considered as a 
lower limit. Additional absorption is ex pected from molecu- 
lar hydrogen. iMauersberger et al.l (1 1995b derived the H2 col- 
umn density in the direction of the nucleus of NGC 253 to 



3.7x10 cm . Taking this value as an upper limit for the col- 
umn density in the disc regions around the nucleus, the column 
densities derived from X-ray spectra are within the limits from 
radio observations. 



4.5. X-ray versus UV morphology 

Fig. [10] shows the X-ray contours from the energy band 0.2- 
0.5 keV overplotted on a two-colour UV image, taken with the 
GALEX observatory (Galaxy Evolution Explorer, a UV space 
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Fig. 10. Two-colour UV image of NGC 253 with NUV (175- 
280 nm) in yellow and FUV (135-175 nm) in blue. The in- 
tensity was rescaled to emphasise the faint, diffuse emission. 
Overplotted are the 0.2-0.5 keV X-ray contours. 

telescope) on 2003 October 13. For the northwestern halo there 
is quite a good agreement between the FUV and X-ray emis- 
sion regions. The FUV emission traces the western horn to a 
distance of ~7.5 kpc above the disc, as well as the broad base 
emission in soft X-rays quite well. In the southeastern halo, 
again, the UV and the soft X-ray emission show the western 
horn structure, where the UV horn extends to about 7 kpc away 
from the disc. However, the FUV horn is slightly offset by 
~700 pc to the northeast with respect to the X-ray horn. Images 
obtained with the Optical Monitor onboard XMM-Newton are 
not sensitive enough to show the extraplanar UV emission. 



Hoopes et al.1 (120051) proposed the following model for the 



origin of the UV emission: Since the UV luminosities are too 
high to to be produced by continuum and line emission from 
photoionised or shock-heated gas, the UV emission could be 
explained by dust in the outflow that scatters the stellar contin- 
uum from the starburst into our line of sight. They also found 
that the UV halo emission, as seen with GALEX, correlates 
with the Ha emission, which could originate from gas that is 
photoionised by UV photons from the starburst. The UV and 
Ha emission would originate in the same cold regions in the 
halo. 

How does the warm gas that is responsible for the Ha and 
UV emission get out into the halo? There are two possibilities: 
either it has already been there from the beginning in the form 
of a cold and maybe clumpy halo component, or it was trans- 
ported by the superwind and galactic fountains from the disc 
out into the halo. 

There are models, where it is possi ble to drag up clouds o f 
cold gas into the halo (e.g. model 3 of IStrickland et al.ll2002l) . 
In a sheet surrounding these clouds, X-ray emission could be 
produced by shocks or in conductive or turbulent mixing inter- 



faces on the cloud surface. This model would also account for 
the non-uniformity of the X-ray emission as seen in the XMM- 
Newton images (Fig. Q] and [3]). However, the model cannot ex- 
plain the displacement of the UV emission in the southeastern 
halo, since in the model the clouds are located within or at the 
inner border of the superwind. 

Could the dust even survive this transport from the 
disc into the halo embe dded in a hot plasma environment? 
Draine & Salpeterl (119791) give the sputtering time for a spher- 
ical dust grain of radius a embedded in a plasma of hydrogen 
with temperatures between 10 6 and 10 9 K and the density «h 
as 

fspu,~10 6 (— I (-^3) V (1) 

\jt/m/ Vcm V 

For wh between 2.5 x 1 0~ 2 cm" 3 in the outflow close to the cen- 
tre JBaueretalJl2007b and 3.2 x 10" 3 cm 4 out in the north- 
western halo, and a grain size of a=0.1 jum, f sput varies be- 
tween 4.0 and 3 1 Myr. So to reach a height above the disc of 
7.5 kpc in less than 31 Myr, an average velocity of at least 
240 km/s is required. This is well compatible with measure- 
ments of outflow velocities in differe nt wavelengths, that ra nge 



from 260 km/s (Na D absorption, iHeckman et al.l uOOOl) to 
about 400-600 km/s (Ha, N II, S II , and Oil emission, lUlrich 
1978c iDemoulin & Burbidgdll970l) . Therefore it is quite pos- 



sible that the dust survives the transport from the disc out into 
the halo. 

Another mod el to explain the UV a nd X-ray morphology 
(e.g. model 5 of IStrickland et al] 120021) requires a thick disc 
component, through which the superwind emerges into the 
halo. On the contact surface between the hot superwind fluid 
and the cold thick disc material we get a heated layer through 
shocks and turbulent mixing where the X-rays are produced, 
surrounded on the outside by a colder layer where the UV emis- 
sion originates. The thick disc component was originally cre- 
ated by lifting material up fro m the disc through the st ar for- 
mation activity (simulations bv lRosen & Bregmanlll995l) . This 
model would easily explain the UV displacement from the X- 
rays, however we would only get a hollow cone with X-ray 
emission. The latter is not what we see in the XMM-Newton 
observations. Though a mix of both models would be able to 
explain the observed morphology. 

For com pleteness, the three other models that 
Strickland et al.l (120020 present are: 



(i) Shocked clouds in the halo (model 2). In this model pre- 
existing neutral clouds are overrun by the superwind, which 
can explain the often asymmetrical patchiness as a result of the 
spatial distribution of the clouds. Ho wever, no HI struc tures 
outside the disc have been found by iPuche et al.l dl99ll) and 



Koribalski et all J1995I) 



(ii) Hot swept-up shell of halo gas (model 4). In this model 
a superbubble bounded by a shock expands into the low-density 
halo medium. The X-ray radiation originates from the shocked 
ambient gas following the morphology of the surface of the 
superbubble. Possible UV/Ha emission could originate from 
dense cool disc gas entrained in the wind. Especially this last 
aspect is not supported by the observations. 
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(iii) Cooling radiation from the wind (model 1). As pre- 
sented bv lStrickland et al.l (120021) this model cannot explain the 
UV or Ha emission, and also does not explain the horn shaped 
morphology of both the X-ray and UV emission. 

The magnetohydrody n amics ISM simulations of 
de Avillez & Breitschwerdtl ( I2005D also shows a clumpy 
halo structure, characterised by turbulent mixing layers, which 
could explain the UV and X-ray filamentary structure. In some 
regions, the magnetic field forms loops surrounded by shells 
which may exhibit enhanced UV emission. 

5. Summary 

We have characterised the diffuse emission in NGC 253. The 
disc extends 13.6 kpc along the major axis and shows emission 
between 0.2 and 10 keV The spectrum could be modelled with 
two thermal plasmas (r co id = 0.1 - 0.3 keV and Th ot = 0.3 - 
0.9 keV) with solar abundances plus a power law component 
and galactic foreground absorption. The power law component 
may indicate an unresolved contribution from X-ray binaries in 
the disc. The total luminosity of the diffuse emission in the disc 
is 2.4xl0 39 erg s -1 (0. 2-10. OkeV). We found clear evidence for 
hot plasma in the disc. The diffuse emission does not originate 
completely from an old stellar population. 

The halo resembles a horn structure which reaches out to 
a projected height of ~9 kpc perpendicular to the disc. The 
halo emission on smaller scales seems not to be uniformly 
distributed, but shows a filamentary structure. The southeast- 
ern halo is softer than the northwestern halo. To model the 
spectra in the halo we needed two thermal plasmas {T co n ~ 
0.1 keV and T^t ~ 0.3 keV) with solar abundances plus 
galactic foreground absorption, which is an indication for the 
non-equilibrium ionisation state of the halo plasma. The to- 
tal luminosity of the diffuse emission is 8.4xl0 38 erg s -1 and 
2.3xl0 38 erg s (0.2-1.5 keV) in the northwestern and south- 
eastern halo, respectively. Densities computed to 3.2 77~ ' 5 x 
10~ 3 crrT 3 and 4.7 77^° 5 x 10~ 3 crrT 3 , with the volume filling 
factor T}. With solar abundances this implies total masses of 
3.3 ?7~ ' 5 x 10 7 M© and 1.8 ?7~ - 5 x 10 7 M for the northwestern 
and southeastern halo, respectively. 

A comparison between X-ray and UV emission showed 
that both originate from the same regions. The UV emission 
is more extended in the southeastern halo, where it seems to 
form a shell around the X-ray emission. 
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Appendix A: EPIC pn images 

We developed an algorithm to create vignetting corrected and 
adaptively smoothed EPIC pn images. This procedure is based 



on the observation itself. Specifically, closed filter observations 
were not used in addition. In the following we will describe the 
algorithm step-by-step. 

The basis for this procedure is a cleaned event file and 
an out-of-time event file. This cleaning included screening 
for high background and also removing bad pixels and bad 
columns. If we want to use more than one observation, all steps 
have to be performed for all observations separately, before the 
products are merged. 

We will here concentrate on the creation of three images, 
which can be combined to a RGB colour image at the end of 
the procedure. In the following all steps are to be done for all 
three energy bands, unless stated otherwise. 

For the desired energy band, we extracted an image from 
the event file. This image was then corrected for out-of-time 
events, via subtraction of an image that was extracted from the 
out-of-time event file and rescaled with the out-of-time event 
fraction. Next, we corrected for the detector background (elec- 
tronic noise, high energy particles) by subtracting the detector 
background surface brightness from the image. This value was 
determined from the corners of the detector which are outside 
of the field of view of the telescope. We here assumed that the 
detector background is uniform across the whole detector (for 
energies above ~7.2 keV this is no longer a good approxima- 
tion, see lFrevberg et alj|2004l) . 

We created vignetting corrected exposure maps and masks 
using the SAS-tasks eexpmap and emask, respectively, which 
will be used to account for differences in the exposure times 
in the images, and to mask the images to regions with an ac- 
ceptable minimum exposure time. Before we could smooth the 
image adaptively, we had to create a template with smoothing 
kernels. This template guarantees that images in different en- 
ergy bands are smoothed with the same kernel size. Therefore 
we added up the images in the different energy bands. With 
the task asmooth, we created the template using the merged 
image, the merged mask, and an exposure map. 

Now we have all the necessary products to smooth the sin- 
gle images. We smoothed the images in the different energy 
bands with the smoothing template, the vignetting corrected 
exposure map, and with the corresponding mask. This step in- 
cluded the vignetting correction via the exposure map, and a 
masking of the image to a region with an acceptable minimum 
exposure time. 

In a final step we used ds9 to create the RGB colour image. 
The resulting images of NGC 253 and a combined RGB colour 
image of the lowest three energy bands is shown in Fig. Q] and 
® 

Appendix B: Background spectra 

The conventional way to create a background spectrum is to 
select a region from the same observation where there is no 
emission from the source. Additionally, the region should be 
close to the source. This way, the spectral background should 
have the same characteristics as the background at the source 
region. In NGC 253, a region which suffices the first criterion 
can be found at the border of the field of view in the southwest- 
ern part of the detector. The second criterion, however, is not 
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satisfied. The background region may show a different detector 
background, and additionally the vignetting is different. 

Since we were interested in determining the characteristics 
of emission with low surface brightness, that extends over a 
large region, where the background is (probably) the dominant 
component, we needed a very accurate estimate of the back- 
ground. Given the very soft nature of the emission, we cannot 
use blank sky observations that were taken in regions of the 
sky where the foreground Nw is different (not to mention other 
uncertainties due to different detector settings, particle radia- 
tion levels, etc.). Here we describe a method to use a local es- 
timate of the sky background that takes properly into account 
vignetting and detector background issues. 

To remove the detector background, we used archival ob- 
servations which were taken in the same mode as the NGC 253 
observations, but where the filter wheel was closed. To avoid 
effects due to changes in the detector settings, or changes of the 
detector performance due to other reasons, we chose the closed 
observations to be as close as possible in time to the NGC 253 
observations. The closed observations we used for observation 
1 and 4 are: revolution 59, obs. id. 0122320701, exposure S003 
(50.5 ks) and revolution 732, obs. id. 0160362801, exposure 
S005 (38.6 ks), respectively. To ensure, that there are as little 
as possible differences between the source observation and the 
closed observation, we removed bad columns and bad pixels 
both in the NGC 253 and closed observation. Additionally, the 
closed observations may have been taken when the spacecraft 
was exposed to a different particle radiation level than the one 
present during observations 1 or 4. The XMM-Newton house 
keeping file contains information on how many CCD rows per 
time unit were rejected due to a possible minimum ionising 
particle (MIP) event, which is a direct estimator of the average 
radiation level. We used these values (see Table lB.il . to rescale 
the count rate of the closed observations. 

We used Out-of-Time spectra from the source and back- 
ground region to correct for contribution from Out-of-Time 
events. When one subtracts a closed observation spectrum from 



Table B.l. Rejected CCD rows due to MIPs per time unit in 
the used observations. 





Obs ID 


Filter 


rejected 


line counter value 


0122320707 
0125960101 


Closed 
Medium 




181.7 
190.2 


0152020101 
0160362801 


Thin 
Closed 




141.4 
120.4 



a Out-of-Time corrected spectrum, one actually removes the 
Out-of-Time spectrum of the detector background twice. This 
is corrected in our method by adding again the Out-of-Time 
spectra of the detector background. 

We corrected the background region spectrum for Out-of- 
Time events and the detector background and applied the vi- 
gnetting correction in each energy bin as a function of off-axis 
angle of the source and background spectrum. This gave us the 
sky background spectrum. 

In all of the above steps, different exposure times and ar- 
eas in the extraction regions have been accounted for. Since 
some of the components in the final background spectrum do 
have low number statistics, we used the conservative approx- 
imati on to Poissonian errors ctn ~ 1 + y0J5 + N (IGehrels 



19861) . To avoid unjustified large errors, we roughly binned the 
spectrum before calculating errors. This binning to a signifi- 
cance of > 3 <x in each bin is performed on the conventionally 
background subtracted source spectrum (no vignetting correc- 
tion and no closed observations) and guarantees that the final 
background subtracted spectrum then has a significance in each 
bin of a least 3 <r. The errors were propagated properly and 
were included in the file with the final background spectrum. 
This spectrum can be used with XSPEC as a background spec- 
trum. 

The whole method can be summarised by Eq. IB.ll with the 
following symbols: 

- B(E) is the counts at energy E in the background spectrum 

- B ODS (E) is the counts in the NGC 253 observation 

- Sdst{E) is the counts from the detector background spec- 
trum in the source region 

- B ODS (E) is the counts in the detector background spectrum 
in the background region 

- 5° b ° T (£) are the counts in the Out-of-Time spectra in the 
source region 

- B^ r (E) are the counts in the Out-of-Time spectra in the 
background region 

- f b s is the exposure time in the NGC 253 observation 

- fdet is the exposure time in the closed observation 

- R b s is the rejected line counter values (see Table IB.U in 
the NGC 253 observation 

- /?det is the rejected line counter values (see Table IB. ll in the 
closed observation 

- As is the area in the source region 

- As is the area in the background region 

- V(E,6s) is the vignetting value in the source region, de- 
pending on the offset angle 9 and the energy E 
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Fig. B.l. Comparison between the background substraction on two examples, (left): the spectrum of region 7, (right): the spec- 
trum of region 14. The top panel shows the new method, as described in this paper, the bottom panel shows the conventional 
method, where the raw background spectrum is used, and a correction for Out-of-Time events has been applied. We only show 
the spectra of observation 4 here, since these have the better statistics. 



- V(E, 9b) is the vignetting value in the background region, 
depending on the offset angle 9 and the energy E 

- f is the fraction of Out-of-Time events in the corresponding 
mode of the observation 

A comparison between this new method and the conven- 
tional method, that does not use the vignetting correction nor 
the closed observations, is shown in Fig. [BJ] for two example 
spectra, both in observations 1 and 4. The single background 
components in the source and background region in observa- 
tion 4 are shown in Fig. IB. 21 All figures show counts integrated 
over the extraction region. The counts in the background region 
were rescaled to the source region area to be able to compare 
them to the source spectrum. Also, the counts in the closed ob- 
servation were rescaled to the exposure time and radiation level 
in the source observation. 

The differences between the new and the conventional 
method in terms of the resulting best fits are the following: 
In the majority of the tested cases, an additional power law 
component with F ~ is required for the fit in the spectrum, 
obtained with the conventional method. The temperatures are 
consistent between both methods, but the resulting flux levels 
in the conventional method are higher. Differences in total flux 
values range between 2% and 22%. The effect between the two 
methods is highest in regions with low surface brightness. Here 
the background dominates and a correct treatment is crucial. 



As an example, the difference in flux in region 7 (low surface 
brightness) is 15% and 22%, for observations 1 and 4, respec- 
tively. Whereas in region 14 (high surface brightness), the dif- 
ferences are 2% and 3%. 



26 



M. Bauer et al.: XMM-Newton observations of the diffuse X-ray emission in the starburst galaxy NGC 253 



N(X ?%>$, R^qinn / 



\CC ?,^3, \i-:::)\cr, 14 




Background region 



Source 

Ou t -of -Time 

Detector 

Out— of— Time Detector 



f%*%HH& 



H 



| !+l " „ 






I' 



WflPJKF 1 ^^ 



Source region 



Source 
Out -or —Time 
Oeiector' 
Out— of— Time Detector 




Background region 



Source 

Out-of-Time 

Detector 

Out— of— Time Detector 




CHrMffl 



Lrsrcy |ksV| 






Lrcrcy |kcV[ 



Fig. B.2. The single components that are part of the total background spectrum compared to the source spectrum, (left): Region 7, 
(right): Region 14, (top): components from the source region, (bottom): components from the background region. The single 
components were corrected for areas, exposure time, and radiation level, with respect to the source spectrum in the source region, 
but no vignetting correction was applied yet. 



